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I. INTRODUCTION. 



Skyrmions in the quantum Hall effect at finite Zeeman coupling. 

A. G. Green, I. I. Kogan and A. M. Tsvelik 
Department of Physics, University of Oxford, 1 Keble Road, Oxford 0X1, UK 

Using a self consistent approximation for the spin distribution of Skyrmions in the quantum 
Hall effect, we obtain an effective action for the Skyrmion coordinates. The energy functional is 
minimised for a periodic distribution of Skyrmions with an underlying hexagonal symetry. We 
calculate the phonon spectrum of this lattice and find that near the classical minimum, breathing 
modes of the Skyrmions are strongly suppressed. The resulting equal sized Skyrmions interact via 
a residual Coulomb potential. Neglecting coupling between phonons and spin waves due to dipole 
and higher order Coulomb interactions, the Skyrmion crystal has a phonon spectrum identical to 
that of an electronic Wigner crystal. The transition to a liquid of equal sized Skyrmions is discussed 
\& , using the theory of dislocation melting and a comparison is made between the predictions of this 

ON ■ theory and the results of a recent experiment. 

0\ 

bJO' 

< 
in 

The theoretical prediction of Skyrmion spin textures in the fractional quantum Hall effect, made by Sondhi et al.B, 
has prompted much theoretical work, supported by mounting experimental evidence, (see Refs. [ Q-Q). 

Here, we consider the low temperature properties of a multi-Skyrmion system. Several parallel approaches have been 
developed to solve this problem. Historically, the first method, used in Ref. employed a field theoretic description 
using a non- linear sigma model with a topological term. A solution for the partition function of topological espitiations 
or Skyrmions in this model, in the absence of Zeeman and Coulomb interactions, was found by Fateev et aln. At zero 
Zeeman energies and for sufficiently high temperature, this solution applies directly to the quantum Hall effectQ and 
describes a meron liquid. When the Zeeman energy is finite, the problem is complicated by the interplay between 
the Coulomb and Zeeman energies. As originallly pointed out by Sondhi et alJj, this fixes the scale for a single 
1 Skyrmion. In the multi-Skyrmion case, the problem is complicated further by the non-local nature of a Skyrmion spin 
distribution. The common belief, based on the analysis of a single Skyrmion, is that we have equal sized Skyrmions 
I ' with some kind of residual interaction. We shall demonstrate that this belief is justified in the physically interesting 
' O ' range of parameters. _ 

An alternative approach is to use a Hartree-Fock approximation!! In particular L. Brey et al. in Ref.[ || suggest 
that the zero temperature ground state of a multi-skyrmion system is a periodic structure rather like a Wigner crystal 
and pick out a face centred square lattice as the most stable of these. Our results, however, indicate that the ground 
state is closer to (or even perhaps, coincides with) a hexagonal structure. 

Recently, a full microscopic.-quantum mechanical solution of the problem with a delta function interaction was 
presented by MacDonald et aln. This approach requires computer diagonahzation of the truncated Hamiltonian. It 
seems that it is rather difficult to apply this approach to systems with large numbers of Skyrmions. 

In this paper we continue to develop the field theoretic description of the problem. In Section II, we introduce the 
energy functional and describe the Skyrmion solution. In Section III, we argue that for small Zeeman coupling the 
Skyrmions are only slightlAz-,distorted by the Coulomb and Zeeman interactions and are well approximated by the 
Belavin-Polyakov solutions^!. We use this solution to investigate the experimentally relevant regime of Skyrmions 
of a small size. In Section V, we show that near the classical energy minimum the Skyrmions all have equal size 
and interact via a residual point-like Coulomb interaction. At low temperatures a Wigner crystal with rectangular 
symetry is formed. We give an analytic expression for this structure and highlight its relationship to the square lattice 
predicted by Brey et al. and the hexagonal lattice expected for an electron Wigner crystal. In section VI, the phonon 
spectrum of the Skyrmion lattice is calculated followed by the spinwave spectrum in Section VII. In Section VIII, 
we discus the dislocation mediated melting of the Skyrmion crystal to form a liquid of equal sized Skyrmions. The 
results of this analysis are compared waith recent experiment. Finally in Section IX, we present a summary of our 
main results. 



II. SKYRMIONS IN THE QUANTUM HALL EFFECT. 
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A. The Energy Functional. 



In ref.[ fjj], the following energy functional is suggested for the ferromagnetic order parameter in the fractional 
quantum Hall effect: 

E Enla Ezeeman ^Coulomb 

= J d 2 x [y |^n| 2 + £2 - g L pBn z ] + ±jd 2 xj d 2 x'q{r)V{r - r')q{r'), (1) 

where, n(r) is an 0(3) order parameter for the spin, p s is the spin stiffness, q the deviation in number density of 
electrons from uniform filling, p is the average density of electrons, v the landau level filling fraction, p is a chemical 
potential for Skyrmions, and V(r) is the Coulomb interaction potential. The filling fraction and average electron 
number density are related by, p = vjlitl , where / is the magnetic length. Aside from the Coulomb interaction this 
is the usual continuous approximation for the energy of a ferromagnet. The Coulomb term arises because, in the 
quantum Hall effect, the electric charge density is proportional to the topological charge density^: 

q = (n.^n x d v ii]) . (2) 
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B. Skyrmion Solution. 

It is useful to parametrise the 0(3) field, n(r), as follows: 



n x + \n v , 2w 



w = - — — , n x + m y = - ; L--l2 , n = , , L , J2 - (3) 



Using this parameterisation of n, the non-linear sigma part of the energy functional is, 



Ad 2 x 



(i + M 2 ) 2 



p s (\dw\ 2 + \dw\ 2 )--^{\dw\ 2 -\8w\ 2 ) , (4) 



where z = x + iy, d — d/dz and d = d/dz. This is minimised on configurations called instantons or Skyrmionsll3, 
which have the form, 

N / - \ 

-w-*n(^|). < 5 > 

i— 1 x ' 

where a,i,bi,h are parameters, aj and hi are called the coordinates of the instanton. According to ref. [ the energy 
of such a configuration is proportional to its topological charge, and does not depend on {h, en, bi}: E = (4np s — p)N. 

III. THE EFFECT OF COULOMB AND ZEEMAN INTERACTIONS. 
A. Distortion of Skyrmions. 

For a particular configuration of Skyrmions, the weak Zeeman and Coulomb interactions cause the spin distribution 
to distort slightly in order to minimise the energy. The energy of these distorted Skyrmions is not degenerate in their 
positions. The original distribution, w — Wo, becomes w = wo + /; Wo(z) is analytic and f(z, z) is some non-analytic 
function which we assume to be small compared with wq. 

Expanding the variation of the non-linear sigma model part of the energy to first order in /, and using Qq = 
ln(l + |wo| 2 ) for ease of notation, we obtain, 

= 8p s [2wdw8w - (1 + \w\ 2 )ddw}/{l + M 2 ) 3 = -8p s B(e- 2Qo df). (6) 

ow 

The total energy is minimised, that is &E 5 nl <* + S(E Coulorn b+E z ^ marl ) _ w j ieil) 
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e- 2Q °df= I dew Q {e)G{w {z),w Q {e)), 



-2Qo; 



df = dew (e)G(w (e) , w (z) ) , 



G(w(z),w (z )) = — e w — . (7) 

8/9 s 

where Q(x) = ln[l + |w| 2 (x)]. The total energy is then, 

E[w] = E nla [w ] + E z+C [w ] + E% a [w ] + Ef +c [ WQ l (8) 

where, 

2Q| 



E% a = 8p s J d'xe-^\df\\ 

E? +c = Jd* X e-Q(wf + wf) S -j^. (9) 

By integrating by parts one may show that E e ^ +C = —2E^ a . From Eq.(|7]) we see that / is of order <?l-B/ p s and that 
the contributions to the total energy in Eq.(||) are of order p s , g^B and (g^B) 2 / p s respectively. 

Therefore, for small Zeeman coupling, p s 3> gLB, we may simply ignore these distortions (therefore from now on we 
shall drop the subscript 0), and minimise on the space of Skyrmion solutions given by Eq.(||). The energy functional 
with these approximations is, 

E = (47r/9 s — p)N + E Zeeman + Ecoulomb, 
Ezeeman = ~9LpB J d 2 x{\ - 2e~ Q ), 

Ecoulomb = \ (U) 2 J d 2 x J d 2 x'AddQ(l)V(x 1 ~ X2 )4dBQ(2) 

= 8 J d 2 Xl J d 2 x 2 VjV( Xl - x 2 )Q(l)aaQ(2). (10) 

In these expressions, Q(l,2) = ^(xi^) and the notation has been relaxed so that w now represents the undistorted 
Skyrmion spin distribution. We will use this convention from now on. 

Taking the functional derivative of Eq.(|l0|) with respect to w, we obtain the extremal condition, 

g L pBe~Q = (^) 2 j d 2 x 2 V\V{^ 1 - x 2 )d<9Q(2). (11) 

We have not been able to solve this equation analytically. It turns out, however, that for the experimentally relevant 
range of parameters we can come up with reasonable approximations. 



IV. APPROXIMATING THE ENERGY FUNCTIONAL AND ITS DERIVATIVES. 



A. The Energy Functional. 

In order to make a useful approximation to the energy functional, we assume that the spin distribution is very 
sharply peaked around the zeros, at. This implies that the topological charge density at the zeros is much greater 
than the average charge density, which is in turn much greater than the charge density at the poles. That is, 
q ai » q » qbi- Later, we will demonstrate the self-consistency of this approximation for the experimental range of 
parameters. 

For spin distributions of this form, we may approximate the Coulomb and Zeeman energies by expanding w around 
its poles and zeros and integrating the corresponding approximate energy densities over patches around these points. 
The leading order behaviour of w near a pole or zero is, 

w(z ~ at) = y/^q a% (z - Oj), 
l/w(z~bi) = y/mi(z-bi). (12) 
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q ai and q bi are the topological charge densities at a, and hi, respectively: 

IQ 12 N N 

nQai = n+k!n 2 H^KH 2 = H 2 n k-^f/nk-M 2 , 



l ' Jll/ "')| 2 w = |aKa i )r 1 | 2 = T ^ fj |6*-M a /ni6.-«il a - (13) 



(1 + 1/M 2 ) 2 1 1 v VJ 1 |b 



The variations of the topological charge density near and hi are, 



nq(z ~ a,) 



(l+7T (Za Jz- ai |2)2 : 



7rg(z ~ 6j) = — ^f' , m2 - (14) 

(1 + Trq bi \ z - °i\ F 

We are now in a position to approximate the Coulomb and Zeeman contributions to the energy functional. Firstly, 
the Zeeman energy is evaluated using the above approximations and integrating over circles of radius d centred on 
each a,i and hi. ird 2 is the area per half Skyrmion and so 2nd 2 is the area per Skyrmion. The radius d is related to 
the average Skyrmion density by, 

q = l/2Trd 2 , q= \v/vo-l\p. (15) 

,/From Eq.([lO|), we find, 

E Zeeman = ^ \J2 ^ ^ ^ + <^/ 2 9) ~ J2 2 ^ + 1bj2q) J . (16) 

q V at bi ) 

Secondly, we evaluate the Coulomb energy. From Eq.([lO|) we find, 

1 x - (ev\ 2 Jnq at f „ „ f 1 1^ (ev) 2 , . 

E Co ulomb= o>^(— / d Xxd X2J— ~, y^JT—, rjyl+o hi r i 1 ') 

2 i VW 6 J (l + \xi\ 2 ) 2 {l + \x 2 \ 2 ) 2 2 e\<n - a,j\ 

dimensionless constant~0(l). 

The first term is due to the local charge density at and the second is an interaction between the charge densities 
near Oj and a,j . The latter has the form of a point particle interaction potential, since all the charge is concentrated 
in small regions of radius 1/ \jirq ai about etj. The first term is dominant because, q ai >> q ~ 1/ — a,j\ ■ Note that 
although the dominant contribution depends only upon the local density of charge, this local density depends on the 
Skyrmion coordinates in a non-local way, through q ai . 
The total energy with these approximations is, 



E= 9£Bj2 [2qj qa% i n (i + q a J 2 q) - 2q/q bi ln(l + q b j2q) + 2ay/q ai /2q) +\Y. T^T" 

% hi 

where, 



<1 



_ ^_ir_£L K ^ , 6}nHn . 1|3 , 



V27r/ e (2tt) 2 eZ#B 

In final expression, we have substituted the typical experimental values of the ratio gLBel/e 2 ~ 0.01 — * 0.0 



(19) 



B. Derivatives of the Energy Functional. 



In order to determine the minimal Skyrmion configuration and its stability, we require the first and second derivatives 
of the energy functional with respect to the Skyrmion positions. These may be approximated by evaluating the 
derivative of the total energy, Eq. (|l~0|) , and expanding the resulting integrals as in the preceeding Section. The 
derivative of the energy with respect to the position of a pole or zero, r i: is 



4 



ff = ± E ln (! + 9aJ2g) - (1 + g afc /2g)- 1 - a^q~j2q) /{a k - rj) 

k 

± 9 -f-Y. ( 2 */«»« ^(1 + ^/ 2 ?) - (! + ^mr 1 ) /(h - r,). (20) 

i 

The plus/minus sign is taken when corresponds to a zero/pole. The second derivatives are given by, 

= ± E (2g/g« fc ln(l + q ah /2q) + (1 + g B */«)(l + 9a fc /2g)- 1 + l/2a^/q a j2qj /(a k - n)(a k - Sj ) 
± 9_pB_Y j ^ (2q/q bk ln(l + q bk /2q) + (1 + q b Jq){l + 9bfc /2g)- 1 ) /(&* - n ){b k - Sj ), (21) 

where the plus sign is taken when both fj and Sj correspond to a pole or zero and the minus sign otherwise. For the 
present we have ignored the contribution from the residual Coulomb interaction. This will be reintroduced later. 

V. THE CLASSICAL MINIMUM. 

The classical minimum is attained when the derivative of the energy is zero with respect to all of the Skyrmion 
coordinates. From Eq. (po|) , this occurs when, to logarithmic accuracy, 

-2 In a N 



q ai =q = 2q ^ fo j , q bi ~ 0, V,. (22) 

As we have already mentioned, our approximation requires qo » q which according to Eq.(|Ti^) is achieved for a 
deviation of v from vq of less than 0.2 (qo/q ~ 5 — > 20 for experimental values of gLBel/e 2 w 0.01 — > 0.04 and 
|^o-l| = 0.2). 

Eq.(p2|) gives 2N constraints on the values of q ai and q bi , and so constrains the positions, a,-,6j. For any particular 
arrangement of ai,we may satisfy the above constraint by fixing The number of degrees of freedom is reduced by a 
factor of two. At the classical minimum the intuitive picture holds. Equal sized Skyrmions with w = z/ro, matched to 
the ferromagnetic background, move independently except for the residual Coulftmb interaction given by the second 
term in Eq.(|17|). In the single Skyrmion sector, this is the result of Sondhi et aln. 

At low enough temperatures Skyrmions interacting in this way are expected to form a classical Wigner crystal. The 
calculation of Bonsall et a/Jl3, performed for an electronic Wigner crystal, suggests a hexagonal lattice as the minimal 
solution. We argue that the same calculation is valid for the crystal of Skyrmions. This brings us into disagreement 
with the results of Ref.[ |). However, this is ?simply? because these authors did not consider a hexagonal lattice 
of the proper type. The simplest periodic arrangement of Skyrmions is described by a bi-periodic analytic (elliptic) 
function, w, with two poles and two zeros per unit cell; b, b' and a, a' respectively. A general property of such an 
elliptic function is that \w(z — a)\ 2 — \w(z — a')\ 2 and \w(z — b)\ 2 = \w(z — b')\ 2 . The dependence of the energy 
functional on |ui| 2 and so n z alone, allows us to treat a and a' , (and also b and 6'), on an equal footing. The Wigner 
lattices have a's and a"s on alternate lattice sites. The hexagonal lattice predicted here, therefore, is not the higher 
energy hexagonal lattice investigated by L. Brey et al, but rather a distortion of their face centred square lattice to a 
face centred rectangular lattice with sides in the ratio 1 : Figures 1 and 2 show the azimuthal components of the 
spin distributions corresponding to the face centred square and face centred rectangular lattice. Analytic expressions 
for these are given by Jacobi elliptic functions; w — h cn(z, l/s/2) for the square lattice and w = h cn(z, e l2jr / 3 ) for the 
hexagonal/rectangular lattice. Since our analysis has neglected dipole and higher moments of the charge distribution, 
the real structure is likely to be a rectangular lattice with the ratio of sides close to 1 : \/3. 

To be sure of this result, we must investigate its stability by evaluating the phonon spectrum. In addition to 
the usual phonon modes of the Wigner crystal there will be modes corresponding to fluctuations in the size of the 
Skyrmions and to spinwaves. 



VI. PHONON SPECTRUM OF THE SKYRME CRYSTAL. 
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A. The Wigner crystal. 

Here we present in condensed form the results of Bonsall et al. for the normal modes of a classical Wigner crystal. 
The residual coulomb interaction is, 



HE-r^. (23) 



and its derivatives with respect to the positions of the Skyrmions are, 

d 2 E 3(ev) 2 a% d 2 E 3(e^) 2 ^ a% 



E 



dciidcij 4e lajd 5 ' daidai 4e 

d 2 E {evf 1 d 2 E {evf 



> j 



15 ' 



E 



ddiddj 4e |ajj| 3 ' dciidai 4e ^-f (a^i 



These are Fourier transformed to give the dynamical matrix, 

7r(ei/) 2 P 



D aa {k) 



2e \k\' 



c G 

where G is a reciprocal lattice vector. To evaluate the normal modes one must first recast these summations into a 
more rapidly convergent form, Rcf. [ [l2| . The result for the hexagonal crystal, written in terms of the longitudinal 
and transverse components, r,s£ {o-l, &t, }i is, 

D rs - y Q T|k| 2 )' L - eVc >' ~ eV c UMba °- (25j 

ao is the lattice spacing and V c is the area of a unit cell. 



B. Contribution of Non-Local Terms. 

For a periodic lattice at the classical minimum, the second derivatives of the nonlocal part of the energy functional 
( ^T| ) are given by, 

dndSj \2irJ 2e ^ {a k - n)(a k - Sj) 1 ' 

Terms like d 2 E/dciidaj are zero for a periodic lattice. Taking the Fourier transform of these derivatives, we obtain 
the elements of the dynamical matrix, 



2eVc gg> V(fc + G)(fc + G') (k + G)G' G(k + G') GG' J ' 



A* = - " E 



GG' 



(k + G)(k + G') {k + G)G l J , 



DM = - " E 



GG' 



<\(k + G)(k + G') {k + G)G'J' 

_^ / e -iG.b e +iG.b \ 



G ^\(k + G) (k + G') 



G 



The summations over the reciprocal lattice vectors, G and G', are only very slowly convergent and must be re-cast 
into a more rapidily convergent form. In order to do this, we use a modification of the Ewald separation method, (see 
refs.[ [nj and [ fL2)). This involves splitting a slowly convergent real lattice sum into summations over the real lattice 
and the reciprocal lattice, both of which are rapidly convergent. Firstly, we write the elements of the dynamical 
matrix in the following way: 

^«=+fe)^=(x:— )(£— J. (3») 

Dss(k) ._ e ^( ? £^)( ? £^). 



Next, using the integral representation, 

1 



Ixl 2 







die"* 1 * 1 , (32) 

and the two dimensional theta-function transformation, 

5^exp{-i|x-a,| 2 +ik.a,} = ^ exp {-|G - k| 2 /4t - iG.x} , (33) 

a; c G 

we consider the lattice sums, 

< k » „ r°° ....... „ m- n ... n . ,, .. 



dte -t|x| +! .k.x (34) 



a;#0 1 M a^O^T G 

a ik.(a(±b) 



^ |a, ± b| 2 -2.7, + l,yj t 



e -|G-k| 2 /4tTi(G+k).b^ (35) 



In these expressions, = q^ 1 is the area of a unit cell and 7 is the Ewald separation parameter. The result is 
independent of the value of 7, which is chosen so that both the lattice and reciprocal lattice summations on the right 
hand side are rapidly convergent. For small wave-vector, k, the G = term dominates; 



e ik.ai 



Yl T^p = V c (- constant - M|k| 2 /47) 



a;#0 



e ik.(a,±b) Trcosrk.b) , .. 10 . \ 

E [a;±b |2 = 7^ (-constant - ln(|k| 2 /4 7 ) . (36) 

We have used the fact that since b lies at the centre of the unit cell, the second summation is independent of its sign. 
Finally, we differentiate with respect to k = k x + ik y to obtain, 

x ^ e ik a ' d ^ e lk a ' _ 2m 1 

^ n ^r~ _ l dk 4- h^" 

ai/0 a, t^O 

^ ai ±b~ dk^\ &l ±h\> v c k ■ 1 ' 

Substituting for these expressions into Eqs.(|28|) to (|3l|), the contribution to the dynamical matrix from the non-local 
interaction is, 
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Dl 





( 








1 


_ e ik.b \ 


Local 










_ e -»k.b 


1 






1 


_ e ik.b 












_ e -ik.b 


1 





o / 



V|k| 2 , 



(38) 



where, 



3q 2 (ev) 



2c 



(39) 



This does not depend upon the lattice type. The interaction is so long range that it is insensitive to the small scale 



distribution of the Skyrmions. 
the derivatives. 



In fact our final result is just that obtained using a continuous Fourier transform of 



C. Phonon Spectrum of the Skyrmion Crystal. 

To detemine the spectrum of the Skyrme crystal, we require the eigen-modes of the total dynamical matrix, 
jjWigner _|_ jjnon-Locai _ -y^ e w jyj concen t r ate on the hexagonal lattice, since this is the minimum energy at our level of 
approximation. In terms of the longitudinal and transverse components, r,s£ {il, o-t, ^t}, the dynamical matrix 
is, 



Dj-a 



£|k| 



T|k| 2 





V 



1 -e 

-j — zk.b 

o l 
l 



k.b 






c ik.b 



-ik.b 



2S/\k\ 



(40) 



The first matrix is given by Eq. (|2q) . The second matrix is obtained from Eq. 
by matrices A T and A respectively, such that, 



by multiplying on the left and right 




/ 1 1 0\ / a\ 

°. ; 01 A\ b \ 



(41) 



-i i 



A- 1 

The problem separates into two independent sectors. In the longitudinal sector, the eigen values and eigen-vectors 
for small wave- vector are, 



CM/2, 

In the transverse sector the eigen-modes are, 

T|k| 2 /2, 



1 

-ik.b 



1 

-ik.b 



and 4<V|k| 2 , 



and 4<5/|k| : 



1 

-ik.b 



1 

,-ik.b 



(42) 



(43) 



The l/|k| 2 modes are extremely stiff, and basically ensure that at finite Zeeman coupling the Skyrmions have a fixed 
size. Note, however, that according to Eq.(|39|) S oc {g^hxg^) 1 ^ and so at zero Zeeman energy we expect substantial 
variation in the size. The remaining modes are the usual Wigner crystal modes. 



D. Quantum Fluctuations of Phonons. 



The above calculation of the derivatives of the energy functional describes the properties of the Skyrmion crystal 
in regimes where the quantum fluctuations may be neglected. We require a full quantum treatment of the phonon 
spectrum. The dynamical part of the action for Skyrmions in the quantum Hall effect is given by the Wess-Zumino 
term, 
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Swz = V 4- I d 2 xdt A[n(r)].?£ = ^ f d 2 xdtdTn.(d t n x d T n) = 2vp [ (fxdtdr ^^ , , (44) 

2 j at 2 j J (l + 

where T = < + ir,dr = d/dT and <9t = d/df. We have extended the function n(i) to a function n(t, r) in the usual 
way, such that t and < r < 1 parameterise the interior of the region bounded by the curve n(i) with n(t, 1) = n(£)EZl. 
This puts the action into a manifestly gauge invarient form. The third expression is obtained by substitution of the 
linear map, Eq.(|^), into this gauge invarient form. 

We may approximate this contribution to the action in the same way as the energy functional, by assuming that 
the spin distribution is sharply peaked around the zeros, ai. Integrating over patches around these points we obtain, 

Swz = 2vp \ d 2 xdtdr (\d T a,i\ 2 - \drai\ 2 ) ^a^x) 
= 2itvf^ I dt dr {Idra,] 2 - |<9 T az| 2 ) 
= -2 W ^/»(a T a i ^a j , s -9 Tfli , 9 ^) 
= -^Jdt(a i , x d t a i , y -a i , y 3 t a i , x ) 

= nvp^ J dtaiedtai, (45) 

i 

where e is the two dimensional anti-symetric tensor, 

«-(;-„'). m 

This is simply the result presented by M. StoneQ for the single Skyrmion and is the Lorcntz force for a particle of unit 
charge moving in a magnetic field 27wp = ev 2 B. Combining this with the approximation to the energy functional 
made above, we obtain an effective action for the phonon modes of the Skyrmion crystal as follows: 



Seff = 5 E / (0^' ~ k ) (iue^Be - D(k)) 5a(«>> k). 



(47) 



Since we are only interested in phonon modes which do not change the size of the Skyrmions, we retain only the 
corresponding elements of D rs . Therefore, D(k) is the two by two matrix £)^* sner /2 gi ven by Eq.(p5|). The equation 
of motion for these fluctuations is 

(iujev 2 Be - D(k)) 8a(uj, k) = 0. (48) 

If we follow Refs.[ |l4|,[l5) and introduce a kinetic term, m*\d t n\ 2 /2, into the Lagrangian, the equation of motion 
becomes, 

(uo 2 m* + iuev 2 Be - D(k)) 5sl(lu, k) = 0. (49) 



Substituting explicit expressions for the matrix -D(k) from Eq.(|25|) we find that the frequencies of the longitudinal 
and transverse modes are, 

^ = 7^|k| 3/2 , 
2uj c m* 

eB . . 

uji^uj c = —-. (50) 
m*c 

This is identical to the result for an electronic Wigner crystal in a magnetic field given in Ref. [ |l3| . 
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VII. SPIN WAVES IN THE SKYRMION CRYSTAL. 



The analysis up to now has been concerned with finding the groundstate spin distribution. We have argued that this 
is a slightly distorted version of the Belavin-Polyakov Skyrmion solution and have minimised the energy functional on 
the space of these solutions. When considering the dynamical properties of the system one must include the possibility 
of spin waves. These have been excluded up to now for a fundamental reason; one cannot write a spin wave in terms of 
distortions of the parameters {ai,bi}. The spinwaves are non-analytic distortions of the BP Skyrmions. For example, 
in the zero charge sector the BP solution is simply the ferromagnetic ground state and does not allow for spin waves. 

Including dynamical terms, the action for Skyrmions in the quantum Hall effect istl 

S = J d 2 xdt^-A[n(r)].^- - J dtE, (51) 

where E is given by Eq.(||) and A[n(r)] is the vector potential of a uni t mo jiopole in spin space. Ignoring the Coulomb 
interaction, this gives the following equation of motion for spin waveal3il3: 

y — =/) s nxV 2 n + pg L nxB. (52) 

These are known as the Lifshitz equations. For small Skyrmions, we expect essentially the same spin wave dispersion 
as for the ferromagnetic groundstate. In fact, on large length scales, the spin distribution looks like a ferromagnet, 
aside from at isolated points where the Skyrmions lie. The spectrum is, 

u{K) = —_{pM 2 +P9lB). (53) 
up 

In addition to the gapped ferromagnetic spin wave, there is an anti-ferromagnetic mode with— a linear spectrum, 
asociated with the anti-ferromagnetic order of the in plane components of the spin distributionlla. In principle one 
may calculate the spectrum of this mode by coarse graining the action for small fluctuations about the Skyrmion 
crystal on a lattice with half the period of the Skyrmion crystal. Staggering the effective action on this lattice, 
taking the continuum limit and integrating out the fast modes will give an effective action for the anti-ferromagnetic 
fluctuations. The solution of this problem is extremely complicated and depends on all the parameters of the system 
in a complicated way through the ground state configuration. In any case, our primary goal is to discus the melting 
of the Skyrmion lattice, a proccess which is unaffected by the presence of spin waves. 

Including the Coulomb interaction will cause a mixing between these spin waves and the Wigner crystal modes 
determined above. It is difficult to estimate this mixing and we will assume it to be negligable. 



VIII. MELTING OF THE SKYRME CRYSTAL. 



There are three steps in our analysis of the melting of the Skyrmion crystal: Firstly, we discuss the regime in 
which quantum fluctuations may be neglected. Next, we describe the various possibilities for defect mediated melting 
of the crystal and compare the predictions with experiment. Finally, we estimate at what temperature significant 
fluctuations in the size of the Skyrmions are to be expected. 



A. Quantum Fluctuations. 



In the usual way, quantum fluctuations of a particular mode may be neglected if the temperature is sufficiently 
high that the first Matsubara frequency, 2ttT, is larger than the energy of the mode at the momentum cut-off. With 
a high momentum cut-off at the first Brillouin zone, |k max | 2 = 7r 2 /a§ = v / 37r 2 g/2, using the phonon dispertion from 
Eq.(po|) and substituting the expressions for £ and T from Eq.(|25[), we obtain the condition, 

T > 0.03^>/^ - 1| 3/2 ~ 5|i//i* - l\ 3/2 K, (54) 

el 

for the temperature above which quantum fluctuations of the phonon mode may be neglected. As indicated, for exper- 
imental values of the paramenters this temperature is only a few tenths of a Kelvin. A similar analysis for spinwaves 
performed in Ref. [ gives the following estimate for the temperature above which these quantum fluctuations can be 
neglected: 
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T>^\u/u -l\, (55) 

which for vq = 1 gives the temperature of order of 10\v/vq — 1\K. If conditions (HJ, ^5|) are satisfied, the Skyrmion 
crystal becomes an entirely classical object. In fact it can be called a crystal only in a loose sense. Since the second 
derivatives of energy are proportional to |k| 2 , the fluctuations of displacements diverge logarithmically : 

2 V C T f d 2 k V C T IT 

where L is the system size. In these circumstances, there is no long range positional order. As usual in two dimensions, 
however, one may have a finite stiffness, T ^ which is associated with a directional long range order. 



B. 2-Dimensional Melting. 

The melting transition in two dimensions is driven by the presence of defects; dislocations and disclinationsEl. First 
attempts at explaining the melting transition treated these two types of defect as independent. The resulting system 
displays two second order Kosterllitz-Thouless type transitions, as first the dislocations and then the disclinations 
unbind with increasing temperature. At the dislocation unbinding transition we pass from the crystal phase, with long 
range directional order and power law decay of positional order, to a liquid crystal phase which has exponential decay 
of potitional order and a power law decay of orientational order, (for an underlying hexagonal symetry, this phase is 
known as a hexatic). The second transition is from the hexatic to an isotropic liquid which has exponential decays 
in both positional and orientational order. This picture, the Kostelitz-Thouless-Nelson-Halperin- Young (KTNHY)cJ 
model, has several successes. Allowing for renormalisations of the elasticities due to the presence of phonons and 
dislocations, the lower temperature transition from the solid to the hexatic phase occurs at the same temperature as 
seen in some experimental systems, (notably the electronic Wigner crystal), and numerical simulations. Unfortunately, 
this is not the full story. Even those melting transitions whose temperatures are predicted correctly show some 
characteristics of being first order. Singularities in specific heat are always too sharp to be adequately described by 
a Kosterlitz-Thouless type transition. j— ■ 

A mechanism by which defect mediated melting in two dimensions could be first order was pointed out by KleinertEll. 
He noted that in fact dislocations and disclinations are not independent. A dislocation can be viewed as a bound pair 
of disclinations and a disclination as a string of dislocations. Although the formation of an isolated disclination has 
an enormous cost in strain energy, this is screened by the presence of dislocations. In fact the dislocation density acts 
like an effective temperature for the disclinations, and virtual fluctuations can lead to a Kosterlitz-Thouless transition 
at which disclinations proliferate. The presence of disclinations also modifies the effective potential of the dislocation 
density, which develops a second minimum at some non-zero density of dislocations. For a certain temperature this 
minimum is at zero energy and the feedback between the dislocations and disclinations drives a first order transition 
directly to the liquid. This occurs at a temperature below the Kostelitz-Thouless transition temperature. This picture 
is still only a caracature of the actual melting process. In fact, dislocations need not pile up on neighbouring sites 
to form a disclinations. If dislocations form a string along next nearest neighbour sites, the resulting configuration is 
identical to a grain boundary and the transition may be driven first order by a proliferation of these boundaries. It is 
sufficient that one allows for the possibility of dislocations forming strings by not choosing a prohibitive form of the 
core energycil. 

The question remains which of these models is applicable to a— particular system. A lattice model of defects 
demonstrating all of these features has been presented by KleinnertEa. Whether the melting proceeds via two second 
order transitions or a single first order transition depends upon the length scale of local rotational stiffness. The 
following partition function describes a lattice with dislocation and disclination defects: 



Z = Il(/ du ^) E exp(-/3£:) 



{riij (x) ,m,i (x)} 

2 



E = - (VjUj + VjUj - l-KUjj) -1 (VjUj - img) + (VjW - 2irnij) (57) 



X,tJ X,l 



The particles are placed at the sites of a simple square lattice with sites labeled, x. Derivatives are all lattice 
derivatives; V/(x) = f(x + i) — f(x),'Vif(x) = f(x) — f(x — i). Ui{x) is the displacement in the i direction of the 
particle at site x, normalised so that when Ui = 2tt, the atom has been displaced by one lattice period, a. The inverse 
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temperature is, (3 = a 2 /j, / (2ir) 2 T . With this normalisation, A and /i = TV c /2 are the usual Lame elastic constants. 
Fure future reference, we define the Poisson ratio as, v = A/ (A + 2/x). Summations over integer mj and even nu 
allow for periodic jumps of atoms to neighbouring sites. This models the presence of dislocations. In the absence of 
dislocations, the first two terms appearing in the energy functional are simply the energy given by linear elasticity 
theory. The presence of disclinations is acounted for by the third term. Disclinations are very singular objects which 
cause large distortions of the lattice. In the vicinity of the disclination linear elasticity theory breaks down, and 
one must include a term which is proportional to the square of the local rotation field, w = (V1U2 — V2Ui)/2. I 
characterises the length scale over which the lattice is stiff to local rotations. 

One final point is that the sum over n,j is unconstrained so that in effect atoms can hop onto neighbouring sites with 
no energy cost. For molecular or electronic crystals, one must modify the theory to allow for the singular repulsion 
between particles. For the Skyrme crystal, this is not necessary since two Skyrmions may easily sit on top of one 
another to form a doubly charged Skyrmion. 

A dual representation of this model in terms of stress gauge fields, may be obtained, by introducing the stress tensor, 
CTy , and the rotational stress tensor, Tj, via a Gaussian transformation and then integrating out the strain fields, Ui(x) 
and w(x). The stress gauge fields are defined in terms of the stress tensors as, <xy = eij\7kAk,n = eijWjh — A{. In 
terms of these, the energy functional takes the form, 



3E = ^{h 



1 



1 1 



2 1 



+ ^7r(Vfc/i - euAif - 2 7 rf(A i & i + h6) 



(58) 



where bi 



and — 'tm, are the integer valued dislocation and disclination densities respectively. 



This partition function may be written purely in terms of the defect fields, bi and 9, by integrating out the stress 
gauge fields, Ai and h. From Eq. ([58]) , the effective partition function in its pure defect form is, 



Z = exp 

{hi*)} 



v.v [ — - v.v 



Vjbj(x) 



(59) 



The character of the phase transitions in this model depend strongly upon the length scale of rotational stiffness. For 
very large I the KTHNY picture holds and there are two continuous transitions. For infinite I and v = 1, the partition 
function reduces to, 



Z = ^2 ex P 



{&«(*)} 



-87r 2 /3^(x)(-V.V)- 



%{x) 



(60) 



which has a Kostelitz-Thouless type phase transition at 4/3 = 2/ir. For very large but finite I, the longitudinal part 
of the field b, is massive and does not contribute to the critical behaviour. The remaining contribution, from the 
transverse component, is given by the first term in the exponent of Eq.(|59|). This may be written, 



Z = exp 

{h(x)} 



-tfW + 1) J2 ^ (x) 6 " V * v J 2 ' Vl >> J (x) 



(61) 



The long range form of the effective potential between the dislocations, is given by, 



<%V.V 



ViV,- 



(V.V) 2 



lim 

5^0 



d 2 q 



(2tt) 2 (q 2 + 5 2 )' 



:(S tj q 2 - qtqj) 



— [Sij \n\x\ - 

47T V J \x\ 2 



(62) 



We have used differentials of the series expansions of the zeroth order Bessel function in order to obtain this results. 
The pairwise interaction energy of the dislocations is, E int = 2tt(3(1 + v) In \x\, and there is an unbinding transition at 
a temperature, (l + i/)/3 = 2/ir, or, T — a 2 fi(l + 1)) / 8ir . This is the universal temperature for the dislocation unbinding 
transition, predicted by KTNHY. In Ref.[ g3|, Morf calculated the renormalisation of the clastic constant, fj,, for the 
electronic Wigner crystal at finite temperature, due to the presence of phonons and dislocations. With a suitable 
rescaling of the parameters, this result may be applied directly to the Skyrmion crystal. The melting temperature, 
allowing for this renormalisation, is 
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r. . „,/„„ _ ir / 2 (^)!_J__ a _ (63 ) 



For small ft ~ a r /I , the effective partition function is, 

Z = J2 cxp ( 2^4tt 2 /3^(9(:z;)(™V.V)~ 1 #(ie) j . (64) 
{00)} \ X j 

This model is equivalent to a Coulomb gas and has a KT transition at 4/3Z 2 /a 2 w 2/-7T, where the disclinations unbind. 

For very small I, there is a single first order phase transition. The constraint V&ft = eki^-i is forced upon the stress 
gauge fields, (see Eq.(§|)). The energy functional then becomes, 



-(VV/i) 2 + 2m (bidj V,/i + Oh) 



(65) 



4/9(1 + v) 

This model shows a single first order phase transition at (1 + v)ft — 0.941. If one studies the same transition in 



the model in which quadratic terms in the original energy functional, Eq. (|57|) , are replaced by cosines, the predicted 
temperature is, (1 + v)(3 = 1.328. Not allowing for renormalisation of the elastic constants at finite temperature, i.e. 
setting /i and v to their zero temperature values, we obtain the following prediction for the melting temperature: 

Allowing for the renormalisation of /i by phonons at finite temperature, [p remains virtually unchanged), this 
temperature is reduced to T m — 0.38|^/^o — ll 1 / 2 . For the cosine model the.-corresponding temperature is, 
T m = 0.34|t//t/o — lp^ 2 without renormalisation and T m = 0.27|^/^q-T l| x / 2 witho. Figure 3 shows a sketch of 
the way in which the renormalisation curve for fj, calculated by Moric3 is used to obtain these corrections to the 
transition temperatures, there will of course be additional renormalisation of the elastic constant which we have not 
allowed for here, due to the presence defects. 

At some intermediate value of I, there is a change from the first order transition to two second order transitions. 
On the basis of numerical simulations to determine the position of this change over, Kleinert has predictedEj that the 
Wigner crystal should undergo a first order transition. The actual value of I for the Wigner crystal lies very close to 
the change over, and so the transition is nearly continuous. We expect that since the Skyrmion crystal shares the 
same phonon spectrum as the Wigner crystal, it too should undergo a first order melting transition. 



C. Comparison with Experiment. 

Recent experimental measurements of the heat capacity of a two dimensional electron system in the quantum Hall 
regime , made by Bayot et alE3 show a marked peak at a temperature of around AOmK with a deviation from total 
filling Vq = 1 of around l^/t'o — 1| ~ 0.2. These authors argue that this contribution to the specific heat comes from 
coupling of the nuclear spins to the Skyrmion spin distribution and that the sharp peak indicates a melting of the 
Skyrmion lattice. The sharpness of the peak suggests a first order transition, which accords with our prediction based 
on the criterion for the local rotational stiffness. The temperature of the first order transition in our lattice defect 
model is a factor of three or so greater than the experimental value at {v/vo — 1| ~ 0.2. There are two points to 
note in reconciling this difference: (i.)\u/uo ~~ 1| — 0.2 is actually outside the regime where quantum fluctuations of 
phonons can be neglected at the melting point. These fluctuations may lead to a reduction in the melting temperature, 
(ii.) Softening of the pairwise interaction between Skyrmions as they begin to overlap means that in regions of large 
distortion such as near a defect, the effective lattice stiffness is reduced and the melting temperature is lowered 
accordingly. 

This latter point may be better seen_by comparing our partition function with that for an N-body system with a 
pairwise interaction potential, $(a; — y)c3. The partition function is, 



Z = 



dx\...dxN 

ivT 



■ exp 



ft 



(67) 



We consider a crystalline state at low temperature, where the particles fluctuate around lattice sites, x, with small 
displacements, u. the partition function may then be written as, 
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a;, 2 a;, 2 



-|^$(x-y + u(x)-u(y)) 



(68) 



Taylor expanding the partition function about the equilibrium configuration, 



z = wXTl 



X.I X.I 



dui{x) I exp 



didj$(x -y)(x- y)k(x - y)ffkUi{x)Viu 3 (y) 



(69) 



For equilibrium configurations including defects, the lattice gradients Vjitj are no longer zero, but equal to the jump 
numbers, riy. Therefore, we make a similar Taylor expansion in powers of WiUj — ria. It is obvious from this final 



form of the partition function and by comparison with Eq. (B7p that a reduction in the curvature of the pair interaction 
potential effectively reduces the elastic constant. In going from the Villain, (quadratic), approximation to the cosine 
approximation for the elastic energy, the non-linearities in the cosine allow for such a reduction in the elastic constant, 
hence the lower melting temperature. 

We conclude that the peak in heat capacity seen by Bayot et al. is consistent with a first order defect mediated 
melting of the Skyrmion lattice to a liquid of equal sized Skyrmions. Up to \v jv^ — 1| ~ 0.05, the melting temperature 
is expected to vary as \v/vq — ll 1 / 2 . Beyond this quantum fluctuations are important and the melting is no longer 
classical. We sketch this behaviour in figure 4. 



D. Fluctuations in the Skyrmion size. 

Now we discuss the criteria for fluctuations in the size of the Skyrmions to be important. The mean square size of 
these fluctuations is 

~W J ~(27rp 64^5 ' (70) 

To estimate when these fluctuations are significant, we use a Lindeman criterion; < u 2 > /V c ~ 1. This is satisfied 
for temperatures greater than 

T > ™ ^f Mn - 1|* (^) " 3 ~ ™0 Ww - II- (71) 

This temperature is always several thousand times the melting temperature and in the experimental range of param- 
eters with \v/vq — 1| ~ 0.2 is around 300.fr. In fact to acheive a realistic temperature for this transition, say of the 
order of IK, would require a reduction in the effective Zeeman coupling by a factor of ~ 100. This takes us outside of 
the regime where our approximations are self-consistent and so we must be cautious about such assertions. However, 
we may be certain that no such fluctuations in the size of Skyrmions occur in the experimental range of parameters 
and that the transition to the meron liquid described in Ref.[ ^| does not occur. 

At temperatures of around 47r / o s significant thermal excitation of Skyrmion anti-Skyrmion pairs will begin to occur, 
and our restriction to distorted BP solutions will be inadequate. 



IX. SUMMARY. 



At moderate Zeeman energies, the distortion of Skyrmions in the fractional quantum Hall effect is small and one 
may describe spin textures using the Bclavin-Polyakov solutions. We show that there is a range of parameters where 
one may self-consistently assume that the charge density of the spin distribution occurs at isolated points. This range 
is spanned by that realised experimentally. The spin distribution of many Skyrmions depends in a very non-local 
way upon the Skyrmion coordinates. However, the interplay of Coulomb and Zeeman energies is such that near the 
classical minimum the behaviour accords with the intuitive picture of equal sized Skyrmions interacting via a residual 
point-like Coulomb interaction. At low temperatures these form a Wigner crystal with hexagonal symetry. When 
one takes note of the azimuthal components of the spn texture, the full hexagonal symetry is no longer apparent; the 
Skyrmions form a face centred rectangular lattice with sides close to the ratio 1 : ^/3. The analytic expression for this 
dstribution is given by w — hcn{z, e 27 ™/ 3 ). In fact, accounting for the higher moments in the Coulomb interaction, the 
actual symetry is expected to be rectangular with sides in a ration somewhere between 1 : y/3 and 1:1. 
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The phonon spectrum of the Skyrmion crystal has four modes. Two of these correspond to the phonon modes of 
the electronic Wigner crystal and two extremely stiff modes correspond to fluctuations in the size of the Skyrmions. 
Quantum fluctuations of these modes may be neglected for filling fractions close enough to vq = 1 and for a sufficiently 
high temperature. We have investigated the defect mediated melting of the Skyrmion crystal in this classical regime. 
A first order transition is predicted which is consistent with recent experiment. Fluctuations of the size of Skyrmions 
and the transition to a meron liquid requires a huge decrease in the effective Zeeman coupling and is not expected in 
the experimental range of parameters. 
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FIG. 1. Square Skyrmion Lattice: plot of (n x ,n y ) derived from w(z) = h cn(z, 1/V2). The dots indicate the zeros of the 
function w(z). The square cell is that chosen by Brey et al. and the parallelepiped is the standard cell. 



t "i ^^^^ 
ft 




FIG. 2. Hexagonal Skyrmion Lattice: plot of (n x 



derived from w(z) = h cn(z,e i27r ^ 3 ). The dots indicate the zeros 



of the function w(z). The solid lines indicate the standard cell and the face centred rectangular cell. 
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H(T)/jl(T=0) 




Cosine model 1 st order 
melting criterion. 

Villain model 1 st order 
melting criterion. 



KTHNY melting criterion. 



renormalisation of \i due to phonons alone. 

renormalisation of u. due to phonons 
and dislocations in the KTHNY model 



FIG. 3. RG flow of elasticity: The solid curve shows the finite temperature renormalisation of the elastic constant \i 
due to the presence of phonons and dislocations. The straight dashed lines show the various melting criteria. The points of 
intersection between the RG curve and the melting criteria give the melting points. Intersections with the horizontal dotted 
line at /i(T)//i(0) give the melting points without renormalisation corrections, see Ref.[23] 




lv/v-11 



.1/2 
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FIG. 4. Predicted melting temperature 

melting temperature and we expect Tm ~ |f /fo ~ 



versus \v/uo — lj: Up to \v/vq — 1| ~ 0.05 the system is classical up to the 
lj 1 / 2 . Beyond here quantum fluctuations may lower the melting temperature. 
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